Localization in Quasi-ID Systems with Random Magnetic Field 



Yakov Rutman, Mario Feingold and Yshai Avishai 

Dept. of Physics, Ben-Gurion University, Beer-Sheva 84 1 05, Israel 



Abstract 



We investigate the localization of electrons hopping on quasi-lD strips in the presence of random magnetic field. In 
the weak-disorder region, by perturbative analytical techniques, we derive scaling laws for the localization length, 
£, of the form £ oc where w is the size of magnetic disorder and the exponent r\ assumes different values in 
the various energy ranges. Moreover, in the neighborhood of the energies where a new channel opens a certain 
rearrangement of the perturbation expansion leads to scaling functions for £. Although the latter are in general 
quantitatively wrong, they correctly reproduce the corresponding rj exponents and the form of the scaling variables 
and are therefore useful for understanding the behavior of £. 

PACS numbers: 71.10,+x, 72.15.Rn 
O^submitted to Phys. Rev. B - June 1995 

On 



00 



Introduction. 



<NThe problem of electron localization in the presence of a 
i^random magnetic field was originally suggested simply as 
^an alternative new type of disorder. 1 Recently however, 
l/~)a realization of the random flux model was proposed by 
^^Kalmeyer and Zhang 2 in order to describe the behavior 
*^of the Hall system in the neighborhood of fractional fill- 
eting fractions with even denominator, v = lp-. While for 

S° ' In 

such f's, the Hall conductance, <r xy , is not quantized, the 

I corresponding longitudinal conductance, <r xx , displays a 
^pronounced minimum which is relatively insensitive to 
^changes of the temperature. This suggests the presence 
Qjof a Hall metallic phase that can be explained by replac- 
ing the entire magnetic flux with flux tubes attached to 
the electrons. If for example, v = then each elec- 
tron will have incorporated two flux quanta, leading to 
a system of composite fermions in zero average magnetic 
field. Moreover, due to the underlying disorder in the 
Hall system, the electron density is fluctuating which 
in turn leads to fluctuations in the associated flux den- 
sity. Correspondingly, a positive magnetoresistance for 
random flux models would account for the observed min- 
imum in <r xx . In turn, the presence of a mobility edge in 
such models is the simplest mechanism leading to posi- 
tive magnetoresistance. In a different context, the quan- 
tum two-dimensional motion of a particle in a random 



static magnetic field is also important for the theory of 
correlated spin system. 3 Sometimes these systems form 
a spin liquid state and such behavior is likely to occur 
in high-T c cuprates. 

At first glance, properties of related models do not in- 
dicate the possibility of a metal-insulator transition in a 
2D random flux model. In particular, in the Anderson 
model with site disorder alone all states are localized and 
the corresponding magnetoresistance is negative. More- 
over, calculations using supersymmetry methods show 
that no extended states are to be found in 2D disordered 
systems with either orthogonal or unitary symmetry and 
only the symplectic case displays a metal-insulator tran- 
sition. On the other hand, in the case of 2D disordered 
systems in a strong constant transverse magnetic field it 
is known that states at the center of the Landau band 
are extended. This is due to an additional topological 
term that appears in the corresponding non-linear sigma 
model Lagrangian, altering the nature of this system and 
allowing for the presence of extended states. In the light 
of recent numerical studies of random flux models that 
found evidence of extended states, 2,4-7 it appears that 
this case is also nonrepresentative of its symmetry class. 

One way of approximating the 2D random magnetic 
field problem is via squares of finite size, L. It was 
found 6,7 that for \E\ < E c ps 3, £,l/L grows as L is in- 
creased while it decreases when \E\ > E c . Such behavior 
indicates that a metal-insulator transition takes place at 
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E c . A different approach to investigating the properties 
of the infinite size 2D systems uses quasi-lD strips of 
length N —> oo and finite width, M . Although all states 
of the quasi-lD system are localized, one expects that 
knowledge of the behavior of the localization length for 
such systems can help to elucidate the above mentioned 
questions. However, numerical results for the random 
flux model on strips 5,6 are to some extent contradictory. 
In Ref. 6 the observed behavior is analogous to that 
obtained for square systems, namely, A(M) = £m/M 
changes from an increasing function to a decreasing one 
at roughly the same value of the energy, E c . On the 
other hand, in Ref. 5 no evidence of a metal-insulator 
transition was observed and A(M) is an increasing func- 
tion for all the energies that were studied. The apparent 
disagreement between the two studies could be a result 
of the fact that, due to numerical limitations, in Ref. 5 
only energies larger than 2.95, | E |> 2.95, were consid- 
ered. A more interesting possibility is that the difference 
between the two studies is related to the short range cor- 
relations that are included in the model of Ref. 6 but 
are absent in that of Ref. 5. 

While all numerical studies were done for large mag- 
netic disorder, extended states are certainly present in 
the limit of vanishing disorder and therefore, one expects 
that this regime is of interest in the search for a metal- 
insulator transition in infinitely wide strips. The advan- 
tage of the approach where weakly disordered quasi-lD 
strips are studied is that a transfer matrix formalism can 
be used for which an appropriate perturbation theory is 
available. 8 Recently, Avishai and Luck 9 have employed 
the perturbation theory of Ref. 8 in order to investi- 
gate the localization of electrons on a ladder network of 
two quantum wires with a random magnetic flux on each 
plaquette. On the other hand, we shall address the cor- 
responding localization problem within the tight-binding 
approximation and for arbitrary M . While in the An- 
derson model disorder is implemented through random 
site energies, magnetic disorder is included via random 
phases in the hopping matrix elements. These phases are 
usually chosen as independent variables that have the 
same probability distribution with some definite width, 
w. The purpose of this work is to obtain a perturbative 
expansion for the localization length, £, in the limit of 
weak magnetic disorder. 

In Section 2, we present our model and derive the 
weak-disorder expansion for the positive Lyapounov ex- 
ponents. The outcome is expressed in several different 
forms corresponding to the various energy ranges. In 
Section 3, these predictions are corroborated by accu- 
rate numerical data for two wires. Finally, in Section 4, 
we discuss the results and some unresolved questions left 
for future work. 



2. The Model 

We study the localization of an electron on quasi-lD 
strips subject to a random transverse magnetic field. 
The model we use consists of a tight-binding Hamil- 
tonian with phase disorder in the hopping matrix ele- 
ments and rigid boundary conditions. The underlying 
2D square lattice has unit lattice constant, o,l = 1, 
length N^oo, 1 < x < N, and width M, 1 < y < M. 
Thus, the lattice sites are located at integer values of 
x and y, n and m, respectively. In order to isolate the 
behavior due to magnetic disorder, in what follows, we 
assume that the site energies, e n>m , vanish, e n rn = 0. 
A random magnetic field, B, with zero average is ap- 
plied parallel to the z-axis, B = B'z. It is chosen such 
that the corresponding vector potential in the Landau 
gauge, A = (0,Bx,0), on each vertical lattice link, 
[n,m] y , between the sites (n,m) and (n,m+ 1), A n>m , 
is an independent random variable. Moreover, all A n>m 
are identically distributed. Although this is clearly not 
the most natural choice, it is among the few models for 
which there is no correlation between consecutive trans- 
fer matrices and accordingly, the perturbative expansion 
of Ref. 8 can be used. Let | n, m > be a complete set of 
orthogonal states associated with the lattice sites. Then 
the Hamiltonian takes the form 

H = | n + 1, m >< m,n \ + \ n — l,m >< m,n | 

n ,m 

_|_ e «a re ,m | „ jm+ 1 >< mj n | 

_l_ e -'««,»-i \n,m—l X m, n | , (1) 

where a n rn = ^j4 n m , and the value of the correspond- 
ing eigenfunction at the (n, m) lattice site, tp n ,m, satisfies 

1pn + l,m + 1pn-l,m + e'""" "* Ipn ,m + l + e~ l " n ' m ~ 1 t/>n ,m- 1 
= Erp nim . (2) 

Eq. (2) can be written in transfer matrix form. Act- 
ing with the 2M x 2M transfer matrix, T n , on the 
2M-vector representing the wave function of two con- 
secutive columns of the strip, t/> n _i m and tp n ,m, where 
1 < m < M , generates the vector corresponding to t^n,™ 
and i/> n _|_i m . The propagation along the strip is therefore 
described by the product 

N 

Q N = n t « • ( 3 ) 

n = l 

and the corresponding localization length is related to 
the Lyapounov exponents, ji , of the infinite product 
Q = limjv^oo Qn- In fact, if Reji > Reji+i, then £ 
is the inverse of the real part of the smallest positive 
Lyapounov exponent. 
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In our model, we can write the random transfer matrix 
T„(a) 



T„(a) 



V n -I 
I O 



(4) 



where I, O and V n are M x M blocks; I is the corre- 
sponding identity matrix, O is the null matrix, and 



I E - e ia ^ o 







E 



V 



(5) 



We assume rigid boundary conditions on the horizontal 
edges of the strip, ip n ,o = tp n ,M+i = 0. Correspondingly, 
the matrix element situated in the upper right (and lower 
left) corner of the V n matrix, vanishes. Moreover, T n 
only depends on the vector potential on the n-th column 
of the strip. Since the a n rn are uncorrelated, this implies 
that the T n for different n are uncorrelated as well, 



T ni T„ 2 = T ni T„ 2 , for n x ^ n 2 



(6) 



where T is the average of the matrix T over the disorder. 
This would not be true if, for example, the Landau gauge 
A = (By, 0, 0) would have been used. In this gauge, T n 
depends on both A n>m and A n -i t m and accordingly, Eq. 
(6) only holds if \n\ — n 2 \ > 1. 

In the site representation, current conservation implies 
that the transfer matrix is constrained to satisfy the re- 
lation 



where 
J = 



J 



-I 

1 o 



(7) 



(8) 



For the case where T n has the structure of Eq. 4, this 
is quivalent to V n being hermitian, V n = Vj, which is 
indeed satisfied by the V n corresponding to our model 
(see Eq. 5). In general, from Eq. 7 one obtains that the 
eigenvalues of complex transfer matrices come in pairs, 
i ^ j, such that, A 8 - = (A^) _1 . Consequently, 
eigenvalues of complex T n that lie on the unit circle are 
at least twofold degenerate. 

In the weak-disorder regime, i.e., when the width of 
the distribution of the random magnetic field is small, 
the Lyapounov exponents can be expanded in a system- 
atic way in terms of the successive moments of this distri- 
bution. This paper presents the perturbative expansion 
of the Lyapounov exponents and gives explicit expres- 
sions for the results up to the second order, following 
closely the approach described in Ref. 8. 



In the case of potential disorder, the corresponding 
transfer matrix is additively modified by the disorder, 
that is, T n = A + ptB n , where /j, is the width of the dis- 
tribution of site energies. For magnetic disorder however, 
this simplifying feature is absent. Nevertheless, one can 
always obtain an analogous separation using the average 
transfer matrix, T n 



(9) 



Accordingly, the zero order approximation to the Lya- 
pounov exponents is given by the eigenvalues of T n , 
namely, jk = logA^. In turn, the A& satisfy the relation 



Afc + X k 1 = E - 2a cos q k 
where 



qk 



id a 



M + 1 



Sk 



(10) 
(11) 



1- 



■ + 0(a 4 ). Moreover, is a permu- 
tation of the integers (1,2,..., M), such that 

I Ai |>|A 2 |>...>|Ajfe |>...>|A M I • (12) 

While for E > 0, Eq. (12) is satisfied if s k = M + 1 - k, 
when E < 0, it implies that s k = k. Since at E = 0, 
| X k |= 1 for all k, the different labeling between the 
positive and negative energy sectors of the band does not 
lead to discontinuities in | A& | at the band center. The 
real part of A^ however, ReA^, flips its sign at E = 0. 
One can see from Eq. (10) that the eigenvalues of T 
depend on a and therefore, the energy where the k-ih 
channel opens, E k , depends on disorder, E k (a). It is 
determined by the requirement that A^(a) = ±1, leading 
to Ek(a) = 2acos</fc ± 2. Therefore, the width of the 
band, A = 2E\, gradually shrinks as the variance of 
the magnetic disorder is increased. To first order in the 
variance, A = 4 + (4 — 2a 2 ) cos M w +1 ■ 

In order to obtain the higher order terms of the per- 
turbation expansion, one must write the random matrix, 
(T n — T n ), in the basis where T n is diagonal. For simplic- 
ity, we assume that the distribution of a n m is symmetric 
and as a consequence, the 0(a 3 ) term vanishes. Since 
the 0(a 4 ) term is rather involved, we only quote here 
the series up to 0(a 2 ) 



p p 

EE 



4a 2 



(M+ l) 2 

T- ■ 



^(A 8 (l)-A 8 (l)-i)(A J (l)-A J (l)-i) 



(13) 



where 
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M 



Tij= sin ^yj ^sin(%)sin(fcgj)[sin((fc - l)qj) 

k = l 

x cos((fc - - sin((fc + l)qj) cos((k + i)^)] . (14) 

It turns out however that the expansion in the moments 
of a is not well behaved for all energies^ E. The most 
serious difficulty originates with the 0(a 4 ) term that we 
ignored in Eq. (13). In Ref. 8, it is shown that this term 
includes sums of ratios of different eigenvalues of T, e.g. 



1_ 

~N 



l<a</3<N 



f)-a 



(15) 



where p < i < M and 1 < j < p. In the limit N —> oo, 
Sij does not converge whenever | A 8 - | = | Xj |= 1 and 
for energies within the band there always exists some 
pair for which that is the case. The only exception 
occurs for the energy range in which there is a single 
open channel, 2a cos + 2 <| E |< 2a cos M w +1 + 2, 

where the 0(a 4 ) term is convergent despite the fact that 
| \ M | = | Xm+i |= I- Another type of divergence occurs 
in the 0(a 2 ) term at the energies at which Afc(l) = ±1 
where 1 < k < p (see Eq. (13)), E k = 2cosg fc ±2. While 
in the case of potential disorder, these divergences oc- 
cur whenever a new channel opens, Ek = Ek, here these 
do not coincide, Ek — Ek = 2(a — l)cosg;% = 0(a 2 ). 
In order to avoid both these divergences, it is enough 
that a strong version of Eq. (12), where all the weak 
inequalities, >, are replaced by strong ones, >, be satis- 
fied. In the following we shall refer to this requirement as 
the nondegeneracy condition. On the other hand, such 
condition excludes most of the energy band which is the 
range of physical interest. Two ways have been suggested 
to obtain information on the Lyapounov exponents using 
the a-expansion in the range where the nondegeneracy 
condition fails. In Ref. 8, the energy is allowed to be 
complex, Ed = E + ie, which for a fixed, nonvanishing 
e insures that the nondegeneracy condition is satisfied 
everywhere except for a finite number of isolated points. 
However, since it is not clear how to extrapolate the re- 
sults down to the real energy axis, this approach, while 
conceptually promising, is merely an uncontrolled ap- 
proximation. On the other hand, in previous work, 11 
we have shown that in the neighborhood of the singu- 
larities of the 0(a 2 ) term, Ek, useful information on 
the localization length can be obtained by ignoring the 
ill-behaved 0(a 4 ) term. In other words, certain prop- 
erties of £ can be extracted from the a-expansion up to 
second order and are not modified by the higher order 
terms. Specifically, for small t = E — Ek and w 2 = a 2 , 
the expansion can be written in terms of new variables, 



Xi = w P'V' , that balance the largeness of the diver- 
gent terms with the smallness of disorder. We refer to 
scaling variables and to the divergent terms as res- 
onant. From the structure of the a-expansion one can 
obtain the /?; and r; exponents. 10,11 Moreover, keeping 
only the resonant terms, the perturbation expansion is 
rearranged such that 



JM = U)' 11 f Vl (x 1 ) + U)' l2 f V2 (x2) + 



(16) 



where r\\ < r/2 < ... < r\i < and f Vt are non-divergent 
scaling functions. Since w is small, we naturally restrict 
our analysis to the first term of Eq. (16). Moreover, we 
find the asymptotic form of f Vl (x\) from the requirement 
that for large x\ it matches up with the resonant part of 
the original perturbation expansion. 

Quantitatively however the prediction of Eq. (16) is 
false, although in some cases the corresponding error is 
extremely small, e.g. for potential disorder with M = 
3 and k = 2. We also noticed that the value of the 
theoretical 0(a 2 ) term turns out to be larger than the 
correct value in the several cases that were studied. One 
is therefore tempted to conjecture that the truncated 
a-expansion leads to a lower bound, £t, f° r the true 
localization length. However, for potential disorder this 
implies that states are extended in the corresponding 
two-dimensional case 12 which is known to be wrong. 

We now consider the scaling approach to Eq. (13) 
for the M > 1 case of the random field model. In the 
neighborhood of the outermost singularity, E = E\ , we 
find that r\\ = 1, f3\ = 2, T\ = 1 and from Eqs. (13) - 
(14) 



fvA x ) - 

M-l 



12^ v 

( T Mj 



(M+ l) 2 
I- t Jm) 



2 cos 7r(1+i) cos 7r(1 ~ i) ) 2 
zcos 2( - M+1 ) cos 2( - M+1 )J 



(1 



(17) 



1 



On the other hand , for E{ with i > 1, r\\ = 4/3, f3\ = 
4/3, T\ = 1. The real part of the corresponding scaling 
function is for x < 



3(M + 1) 2 V^W (1 - 2 cos cos $=ft 



7r(l-Q 



1 



For x > however, all the resonant terms in Eq. (13) 
are imaginary. 



(18) 
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3. Comparison with numerical simulations. 



In order to compare the theoretical predictions with the 
results of numerical simulation, we consider a strip of 
width M = 2 and a uniform distribution for the vector 
potentials, a n rn , 



P (^n ,m) 



u 



for 
for 



|< w 1 
|> w l 



(19) 



where w\ = V3w. The spectrum of the ordered system 
consists of three domains, namely, I. | E |> 3, II. 1 <| 
E |< 3 and III. | E |< 1. Their edges are defined by 
the A2(l) = ±1 and Ai(l) = ±1 conditions, respectively 
and the localization length has different forms in each 
of these regions. Since the behavior of the localization 
length is invariant under E —> —E, in what follows, we 
only consider the positive energy range, E > 0. 

In the region 1 < E < 3, one of the A-eigenvalues is 
real, Ai, and the other, A2, is imaginary. In other words, 
only one transversal channel is open. In order to obtain 
72, and hence the localization length, we have used the 
perturbative result of Eq. (13). Since A2 is imaginary, 
neither the term of order zero nor the term proportional 
to w 2 contribute to 72. As a consequence, only fourth- 
order terms contribute to 72 and therefore, the inverse 
localization length vanishes as w 4 . Indeed, in Fig. 1 
the inverse localization length at E = 2.5 is shown to 
grow as wf with k = 4.00 ± 0.03, in agreement with our 
prediction. 



10" 



10" 



10 = 



1Cf 




0.1 



1.0 



Fig. 1. Numerical values of the inverse localization 
length, as a function of the strength of disorder, w\, 

at E = 2.5 (error bars). The straight line is the best fit of 
the form f" 1 = Crof, where k = 4.00 ± 0.03. 
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-0.2 



0.2 



0.6 



Fig. 2. The real part of the scaling function fi(x) de- 
scribing the behavior of smallest Lyapounov exponent 72 near 
external band edge, E = 3. The data correspond to rectan- 
gular distributions of the random values, with various values 
of the width, Wl = 0.1 (A), 0.15 (x) and 0.2 (o). The con- 
tinuous line represents the scaling curve given by Eq. (17). 
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0.00 
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-1.5 



-0.5 



Fig. 3. Same as in Fig. 2 only near internal band edge, 
E = 1. The continuous line represents the scaling curve given 
by Eq. (18). 

For energy E < 1, both eigenvalues, Ai and A2, are 
imaginary, i.e. both transversal channels are open. Un- 
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fortunately, the general weak-disorder perturbative re- 
sult of Eq. (13) does not yield directly any quantitative 
information about the Lyapounov exponents in this re- 
gion. While for large enough e the perturbation theory 
is valid, such addition to the energy changes the physics 
of the problem; the perturbation theory prediction now 
refers to the £ corresponding to a physically irrelevant 
model with complex energy that significantly differs from 
the one with real energy. On the other hand, if the min- 
imal e required for regularizing the perturbation expan- 
sion, e c , is not too large, one could hope that £(e c ) is a 
good estimate of £(0). Unfortunately, numerical simula- 
tion shows that this is not the case. Consequently, both 
Lyapounov exponents vanish like w 2 , but one cannot ob- 
tain explicit expressions for the amplitude functions. 

We now look at the behavior of the second Lyapounov 
exponent, 72, near the external band edge of the spec- 
trum of the unperturbed network, that is, for E —> 3. 
For x —> +00, that is, deep outside the band, the result 
for the scaling law can be written in the form 



72 = w 



(y/x- 



where 



E-3 



Conversely, for x 



(20) 



-00, i.e., deep in 



the region with a single open channel the second Lya- 
pounov exponent vanishes as the fourth power of the 
strength of disorder. Fig. 2 shows the real part of the 
scaling function /1 = obtained from data correspond- 
ing to a narrow distribution of vector potentials. For 
x > 0, scaling is observed to hold extremely well and 
so does Eq. (20) for large enough values of x. On the 
other hand, for negative x, scaling holds only in a small 
interval, up to x ps —0.05. 

For E —> 1, the spectrum of the network exhibits an 
internal band edge, which demarcates the two-channel 
region from the one-channel region. The perturbative 
weak-disorder expansion of Eq. (13) for the second Lya- 
pounov exponent is again singular. In analogy with the 
previous case, we are led to hypothesize the scaling form 



72 



„ 4 /3. 



1 



for 



(21) 



where x = i s the corresponding scaling variable. 

For x > however, all the resonant terms in Eq. (13) are 
imaginary and one cannot obtain information about this 
part of the scaling function from the 0(a 2 )-terms. Fig. 
3 shows the real part of the appropriate scaling function, 
/ 4 / 3 = J^ /3 . While scaling holds here as well as it does in 
Fig. 2, unlike in Fig. 2, the theoretical scaling function 
is significantly larger than the numerically obtained one. 
This is a consequence of the failure of the a-expansion 
in this range. 



4. Conclusions 

We have investigated the localization of electrons in the 
framework of a tight-binding Hamiltonian on quasi-lD 
strips with a random magnetic field. The magnetic vec- 
tor potentials are assumed to be independent and drawn 
from a common even distribution. As in any model of 
ID disordered wires, all eigenstates are exponentially lo- 
calized. We have therefore focussed our attention on 
the weak-disorder regime, w <C 1, where the localization 
length, £ = 1 /Re7M , is much larger than the lattice spac- 
ing, £ d£. In this regime, a systematic perturbative 
expansion for the smallest positive Lyapounov exponent 
was derived. Moreover, we have checked the outcome of 
our analytical approach against data obtained by means 
of numerical simulations for M = 2. 

The most important result is that magnetic disorder 
modifies the divergence laws of the localization length 
£ in the weak-disorder regime (w — ► 0), in the various 
energy domains, with respect to those observed in the 
case of potential disorder. If we go continuously from 
outside the band toward its interior, we encounter the 
following sequence of exponents for the divergence law 
of the localization length, £ oc l/w rl , 



r, = 0,1, 4, 4/3, 2, 4/3, 2, 



(22) 



This sequence is more diverse than the one observed for 
potential disorder 11 , namely 



r? = 0,2/3,2,4/3,2,4/3,2, 



(23) 



We suggest that it should be interesting to use the 
approach described in this work to study the combined 
effects of potential and magnetic disorder. 13 Moreover, 
it would be useful to generalize the theory of Ref. 8 
for the case where transfer matrices corresponding to 
neighboring strip columns are not uncorrelated and find 
the influence of such correlations on the ry-exponents. 

However, the most important task remains finding 
the way to properly regularize the perturbation expan- 
sion for infinite products of transfer matrices such that 
it would become quantitatively correct throughout the 
band. This problem is reminiscent of analogous diffi- 
culties arising in the Hamiltonian and Green function 
perturbation theories and one expects that some of the 
methods developed there could be applied to the case of 
infinite products as well. 
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